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We apply the stochastic variational method to the action of the ideal fluid and 
showed that the Navier- Stokes equation is derived. In this variational method, the 
effect of dissipation is realized as the direct consequence of the fluctuation dissipation 
theorem. Differently from the previous works [5,6, we parameterize the Lagrangian 
of SVM in more general form. The form of the obtained equation is not modified 
but the definition of the transport coefficients are changed. We further discuss 
the formulation of SVM using the Hamiltonian and show that the variation of the 
Hamiltonian gives the same result as the case of the Lagrangian. 
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1. Introduction : origin of stochasticity in variational approach 

The variational principle is one of the important guiding principles in physics. The 
usual variational principle, however, is not applicable when irreversible dynamics is 
present. It is because dissipation is caused by energy exchange processes between macro- 
scopic and microscopic motions whose information is not included in Lagrangian. For 
example, the Euler equation can be derived from the Lagrangian £(£,p m ,v) which is a 
function only of hydrodynamic variables: the energy density e, the mass density p m and 
the fluid velocity v. On the other hand, the viscosity of fluids is induced by the collisions 
between the constituent molecules of fluids. Thus, to introduce the effect of viscosity in 
the framework of the variational method, a new term, which represents the microscopic 
degree of freedom (molecules), is added to the original Lagrangian. This approach is 
known as Rayleigh's dissipation function method, but it is not easy to specify uniquely 
this additional term. 

In the stochastic variational method (SVM) [U[2], this microscopic effect is taken into 
account through a different mechanism. Note that a fluid is approximately interpreted as 
the ensemble of fluid elements inside of which there are many molecules and the thermal 
equilibrium is achieved. The fluid elements move following the flow of the fluid, and the 
collisions among them are the origin of viscosity. Let us pay attention to the motion of 
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Fig. 1: The collisions between fluid elements are interpreted as a kind of Brownian motion 
of a fluid element under the noise of molecules. 

one fluid element as is shown in Fig. [TJ Then, because of the fact that each fluid element 
is consisted from many molecules, the interaction between this and other fluid elements 
(shown on the left panel) can be regarded as the collisions between this fluid element 
and molecules (shown on the right panel). This picture is very similar to the Brownian 
motion. 

Motivated by this fact, the definition of velocity of the fluid element is given by a 
Langevin equation in SVM. Then the information of the microscopic degrees of freedom 
which is not included in the Lagrangian are introduced through the behavior of the noise. 
In this approach, the fluctuation-dissipation theorem is naturally satisfied and the non- 
dissipative dynamics (the Euler equation in the case of hydrodynamics) is reproduced in 
the vanishing limit of the noise. 

So far, SVM has been exclusively applied to re- formulate the Schrodinger equation 
from the viewpoint of the variational method [3]. However, because of this consistency 
between SVM and the fluctuation-dissipation theorem, it is natural to expect that SVM 
is still powerful method to discuss dissipative phenomena. 

For the application of SVM to derive the incompressible Navier- Stokes equation has 
been studied [4]. Recently, it was shown that even the compressible Navier- Stokes equa- 
tion can be derived [5] [6]. In this paper, we generalize the previous argument. We 
parameterize the stochastic Lagrangian of SVM in more general form and introduce a 
Hamiltonian in terms of the Lagrange transform. 

This paper is organized as follows. In Sec. 2, we develop SVM in the Hamiltonian 
formulation and apply it to the derivation of the Navier-Stokes equation. In Sec. 3, we 
discuss that the fluctuation-dissipation theorem is naturally realized in SVM. Section 4 
is devoted to concluding remarks. 
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2. Application to hydrodynamics 

In the following, we develop the formulation of SVM and derive the Navier- Stokes 
equation. The technical details of the derivation are mainly following Refs. |5j[6]. 

2.1. Definition of velocity 

In the classical variational method, the velocity of a particle is a continuous function 
of time. On the other hand, as was mentioned in the previous section, the motion of 
the fluid element fluctuates because of the random collisions with molecules, and the 
trajectory is determined by the following stochastic differential equation (SDE), 



dr(t) = u(r(t),i)dt + V2udW(t). (dt > 0), (1) 

where v is the strength of the noise and 

EiWj^-Wjiu)} = 0, (2) 

EliWjW-WjiuWWkW-Wkiu))] = V k \s-u\, (3) 

EKWj^-WjiuVWkit)} = 0, (4) 



In SVM, the effect of the microscopic degrees of freedom which is not included in the 
Lagrangian is taken into account through the noise dW(t). 

The trajectory r(t) is not smooth and hence we cannot define the time derivative at t 
uniquely, as is shown in Fig. [2] Consequently, we further introduce another SDE which 
describes the time reversed process of Eq. (pQ), 

dr(t) = u(r(t), t)dt + V2udW(i), (dt < 0) (5) 

The noise dW satisfies the same correlation properties as dW and there is no correlation 
between dW and dW. The equation must describe the time reversed process of Eq. (pQ). 
To satisfy this condition, the Fokker-Planck equations derived from Eqs. (pQ) and (j5j) 
should be same. Then u is related to u, as 

u = u + 2zA71np, (6) 

where p is the particle density which is given by the solution of the Fokker-Planck equa- 
tion. 

Correspondingly to the two SDEs, we introduce the two different definitions of the 
time derivative [7]: one is the mean forward derivative, 

Dr = u, (7) 

and the other the mean backward derivative, 

Dr = u. (8) 

The time derivative terms which appear in the kinetic term of the Lagrangian are replaced 
by Dr and/or Dr. 
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Fig. 2: The ambiguity for the definition of the time derivative of r. 



2.2. Definition of action 

As is discussed in Ref. [5], the classical action which leads to the Euler equation is 
given by 



where t a and t& are initial and final times. Note that the action is expressed in the 
Lagrange coordinate and p™ denotes the initial mass density. The relation between the 
particle density p and the mass density p m is given by p m / p = m, where m is the mass of 
the constituent molecules of the fluid. When we apply the classical variational method 
to this action, we obtain the Euler equation as is discussed in Ref. [5|. 

As was emphasized in the introduction, the effect of dissipation can be taken into 
account through the noise in SVM and we still apply the same action for the variation. 
However, because of the two different definitions of the time derivatives, D and D, the 
stochastic representation of this action is not unique. Note that the dissipative (viscous) 
terms should disappear and the Euler equation must be reproduced in the vanishing noise 
limit. That is, the stochastic action should be reduced to the classical action in this limit. 
Then the most general expression of the stochastic Lagrangian density is given by 




(9) 
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where the Lagrangian density in the Lagrangian coordinate is 



L = 
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(11) 

Here ai and «2 is arbitrary constants. In this expression, we assume that the terms 
higher order than (dr) 2 and (Dr) 2 do not appear. One can easily check that, in the 
vanishing noise limit, v — » 0, the difference between Dr and Dr disappears, 

lim Dr lim Dr = dr/dt, (12) 



and hence, independently of the values of oi\ and «2, the above action ([10]) always 
reproduces the classical action. 

Note that this action is more general than that used in Refs. [5 , which is reproduced 
by setting a\ = 1/2 and (\2 = 1/2. 

2.3. Various formulas 

The most part of the variational procedures in SVM is same as those in the classi- 
cal variational method. However we cannot use the Taylor expansion and the partial 
integration formula. These are modified in SVM as follows. 

1. Ito formula [8] 

For an arbitrary function F of the stochastic variable r satisfying Eq. (pQ), there is 
the following expansion, 

dF(r, t)=ftF + u-VF + vV 2 F + 0(dt 3/2 ). (13) 

2. stochastic partial integration formula [9] 

For the stochastic variables X and Y, the partial integration formula is extended 
as follows, 



t b 

dtE[DX(t) • Y(t) + X(t) • DY(t)] = E[X(t b )Y(t b ) - X(t a )Y(t a )}. (14) 



2.4. Stochastic variation 

In Refs. [5] [6], the stochastic variation was discussed based on the Lagrangian. Here, 
we develop the argument using the Hamiltonian which is introduced by the Legendre 
transform of the Lagrangian. 

From Eq. ([TT]) . the conjugate variables are introduced as 



PO 

2 
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where the initial particle distribution is given by po = p™/m. From the Legendre trans- 
form, the Hamiltonian corresponding to this Lagrangian is given by 

«(r )P) p) = ^ (pDr + pDr) -C 

p ° r G + a2 ) { G - ai ) p2 + G + a } - G - a2 ) pp 



8m 2 M 



(17) 



where M = a 2 /2 — afa^ — 1/4. In this derivation, we used that M/0. 
Then the action which is a function of r, p and p is given by 



Ar,p,p) 



dt / d 3 KE 



Po- 



pDr + pDr 



n 



(18) 



Note that this action depends on three variables, r, p and p. Differently from SVM in 
Refs. [5j[6], we have to consider the variations for these three variables. Then we obtain 



Dt 
Dt 



(l/2 + a 2 )(l/2-ai) 1/2 - a 2 



Dp + Dp 



2mM 
(l/2 + a 2 )(l/2 + ai) 
2mM 

m- =0. 

St p m 



4mM 
1/2 - a 2 
4mM 



P = 0, 
p = 0, 



By eliminating p and p from the last equation, we obtain 



0L2 



{a 



a\ I Du 



ai I Du 



a 2 (Du + Du) 



5_j_ 
5r p m 



(19) 
(20) 
(21) 

, (22) 



Here we used Eqs. (pQ) and (j5j) to calculate the forward and backward mean derivatives. 

In general, the energy density e depends on the specific entropy and hence, to calculate 
the variation, we have to need certain model for the variation of the entropy. See Refs. 
[5j El for detail. 

As a result, the equation corresponds to the Navier-Stokes equation is obtained by 
the stochastic variation, 

!«2 



P m (d t + v m • V) w % m + dj((P - CV • v m )5" - V e%) - 4a 2 is 2 p m d z (^ d 2 ^T) = 0, (23) 



where 77 = ai(l + 2a 2 )vp m , e™ = djW l m + diV J m 
is given by the mean of u and u, 



u + u 



v m = 



V • v m and the fluid velocity v m 

(24) 
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In this derivation, we used the model for the variation of entropy used in Ref. [5] and 
then the bulk viscous coefficient £ is introduced. This result is the generalization of that 
of Ref. [5j [6]. When we use ql\ = a 2 = 0, the effects of the noise cancel each other and 
Eq. ([23]) is reduced to the Euler equation by ignoring the effect from the variation of the 
entropy, that is, £ = 0. When we choose ql\ = 0, a 2 = 1/2 and appropriate v and energy 
density e, we obtain the Gross-Pitaevskii equation. See Ref. [6] for details. 

When we set a 2 — and abrbitrary finite a±, this equation is completely equivalent 
to the Navier- Stokes equation. However, in the framework of SVM, the effect of the noise 
induces not only viscous terms but also higher order correction term to the Navier- Stokes 
equation. This correction is not only of second order for the magnitude of fluctuations 
(noise) but also of third order for the spatial derivative V. The Navier- Stokes equation 
does not contain such a term since, by construction, only the first and second order 
spatial derivative terms are maintained. In the case of rarefied gases, this corresponds to 
the first order truncation in the derivative (Chapman-Enskog) expansion of one-particle 
distribution functions. However, there is a possibility that this correction term becomes 
important in the discussion of very viscous fluids and the behaviors of the shock wave. 
Interestingly, this term may be identified with the surface effect of the Thomas-Fermi 
approximation used in the nuclear physics. 

In SVM, the behavior of the mass density p m is automatically determined by the 
Fokker-Plank equation corresponding to Eq. (pQ) (and (|5j)). By using v m , this is expressed 
as 

ftp m + V(p m v m ) = 0. (25) 
This is the conservation equation of the mass density. 

3. Fluctuation-Dissipation theorem 

As was mentioned in Ref. [5] , the noise introduced in Eq. (pQ) is related directly to the 
transport coefficients. For example, we can show that the kinetic viscosity r)j p m satisfies 
the Einstein relation, 

n 1 f°° 

= - ai (l + 2a 2 ) / dtE[Sv(t) • 5v(0)], (26) 

P o Jo 

where 5v = dr/dt — u =\/2vdW / dt. Note that the result in Ref. [5 j is reproduced by 
using ai = 1/2 and a 2 = 1/2. This is nothing but the realization of the fluctuation- 
dissipation theorem and it appears as a natural consequence of SVM. We thus conclude 
that SVM possesses not only the well-defined mathematical structure but also a reason- 
able mechanism of dissipation. 

4. Concluding remarks 

We applied the stochastic variational method to the action of the ideal fluid and 
showed that the Navier- Stokes equation is naturally derived. Differently from the pre- 
vious works [5j [6], we parameterized the Lagrangian of SVM in more general form. 
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The structure of the obtained equation is not modified but the definition of the trans- 
port coefficients were changed. We further discussed the formulation of SVM using the 
Hamiltonian and showed that the same result can be obtained from the variation of the 
Hamiltonian. 

The stochastic Hamiltonian is discussed also in Ref. [10 j, but only one conjugate 
momentum is introduced there. On the other hand, two conjugate momenta are necessary 
in our approach. The existence of the two different momenta becomes important to 
discuss the uncertainty relation in hydrodynamics [TT] . 

What is the importance of the formulation of hydrodynamics in the framework of the 
variational method? If this approach can be established as another principle of variation, 
we can derive various dissipative equations in more systematic manner. As a matter of 
fact, the variational method with the Rayleigh dissipation function has been used to 
derive the evolution equations in soft matter physics [12]. And it may become possible 
to discuss the general properties of different dissipative phenomena from the view point 
of the symmetry of Lagrangian. As for the Noether theorem in SVM, see Ref. [T3] . 

As a matter of fact, in SVM, the effect of dissipation is induced by the effect of 
fluctuations and the fluctuation-dissipation theorem is automatically satisfied. This will 
be an evidence that SVM accounts for the essential mechanism of dissipation, and indi- 
cates that SVM can be applicable not only to hydrodynamics but also to more general 
dissipative phenomena where the fluctuation-dissipation theorem is satisfied. 

As shown, the present result of SVM specifies the form of the higher order correction 
to the Navier-Stokes equation. If SVM is a reliable approach as was discussed above, 
this higher order term neglected in the Navier-Stokes equation should be considered 
seriously. In fact, as is discussed in Ref. j6|, this higher order term plays a crucial role in 
diffusion processes. The application of SVM to a diffusion process leads to a generalized 
form of the diffusion equation which contains memory effects [14], and the higher order 
term guarantees that this generalized equation reproduce Fick's law in a certain limit. 
Moreover, this term becomes significant when the inhomogeneity of the mass density 
increases. Thus the behavior of the shock wave will be affected by this term. 
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